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It has been shown that interband mixing of gradients of two order parameters (drag effect) in an 
isotropic bulk two-band superconductor plays important role - such a quantity of the intergradients 
coupling exists that the two-band superconductor is characterized with a single coherence length 
and a single Ginzburg-Landau (GL) parameter. Other quantities or neglecting of the drag effect 
lead to existence of two coherence lengths and dynamical instability due to violation of the phase 
relations between the order parameters. Thus so-called type-1.5 superconductors are impossible. An 
approximate method for solving of set of GL equations for a multi-band superconductor has been 
developed: using the result about the drag effect it has been shown that the free-energy functional 
for a multi-band superconductor can be reduced to the GL functional for an effective single-band 
superconductor. 
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I. INTRODUCTION 


Two-band superconductors are a specific class of superconductors essentially differing in their properties from single¬ 
band superconductors. Their typical representatives are magnesium diboride MgB2, strontium ruthenate Sr2Ru04, 
nonmagnetic borocarbides LuNi2B2C, YNi2B2C and ferropnictides. In this article we will consider only isotropic bulk 
(polycrystalline) s-wave superconductors. One of the main feature of these materials is the presence of two energy 
gaps Ai and A2 which, however, vanishes at the same temperature (FigQ. According to microscopic theory 
m presence of the two gaps is explained by the fact that in each band i an own coupling constant ga exists - 
the intraband constant. In the same time, the interband coupling constant gij exists too, which, on the one hand, 
enhances pairing of electrons, on the other hand, leads to the single critical temperature Tc. BCS gap equations for 
a two-band superconductor are [TH3]: 

giiAi tanh(F;i,fc/2fcBT) gi2A2 tanh(F;2,fc _ 

, 2£'2,fc 

k 

_ 522A2 tanh(i?2,fe/2fcgT) 512A1 tanh(£'i_fc/2fcgT) 

k ’ k 


k 



where is the quasiparticle’s energy in a band i. Unlike single-band BCS theory a superconducting state can exist 
both attractive interband coupling constant gi2 > 0 and repulsive gi2 < 0, moreover the gaps are nonzero if even the 
intraband couplings are absent gn = g22 = 0. In the case of the attractive interband interaction the gaps have the 
same phases on both Fermi surfaces, while for the repulsive interaction the phases will be opposite. Thus the phase 
difference of the order parameters lAije^^b |A2|e®'^^ are: 

cos((/?i - (P2) = l if gi2 > 0 .2) 

cos(v?i - (P2) = -1 if gi2 < 0 

For example, in absence of magnetic field we can suppose Ai > 0, then we will have A2 > 0 for gi2 > 0 and A2 < 0 
for gi2 < 0. From Eq. 0 we can see the important property of a two-band superconductor: if we violate the phase 
relation ([^ then suppression of the energy gaps Ai, A2 will take place (extremely strong suppression if the intraband 
couplings are absent gn = g22 = 0). For the suppression of the order parameters the violation of the phase-locked 
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states (p = (/?! — (^2 = 0 or TT must be macroscopic, constant in time and not small (for example, when there are 
two different coherence lengths in a system with two gaps), unlike Leggett’s mode, which is collective mode of small 
fluctuations of relative phase ip{v, t) and behaves like the Anderson plasmons in Josephson junctions [51 |B]. In addition, 
we assume that current in a two-band superconductor is less than some a critical current J < Jt over which interband 
phase breakdown occurs, resulting in spontaneous phase solitons in ip{r,t) ( 3 , which is nonequilibrium state. In 
phenomenological theory the coupling between the bands is represented by Josephson-like coupling term: 

+ ( 3 ) 

in a free energy functional, where 'Ll and '^2 are order parameters for band 1 and 2 accordingly. 

Currently there are two opinions about the properties of two-band superconductors: 

1 ) In papers |8l|9j it has manifested about a new type of superconductivity in MgB2 - a novel ’’type- 1.5 supercon¬ 

ductor” , contrary to type-I and type-II superconductors. In papers unHH] a two-band superconductor was studied, 
where they considered GL parameters Ki = X/^i {i = 1 , 2 ) in two different regimes to produce type-I (ki < l/v^) and 
type-II (k 2 > l/'/S) materials, that corresponds to different coherence lengths and ^2 = Ja2(o) ■ That is 

each correlation length is sorted with a corresponding band, where Fermi velocities vpi,VF2 and energy gaps Ai, A2 
are different. Their prediction leads to what they call a ”semi-Meissner state”. Instead of homogeneous distribution, 
the vortexes form aperiodic clusters or vortexless Meissner domains, arising out of short range repulsion and long 
range attraction between vortexes. 

2 ) However, in review [ 13 ] an opposite opinion has been suggested in respect of existence of the type- 1.5 supercon¬ 
ductivity in two-band superconductors. It was shown that for the real superconductor MgB2 which possesses a single 
transition temperature, the assumption of two independent order parameters with separate penetration depths and 
separate coherence lengths is unphysical. In particular, in the above-mentioned works [ii numerical estimates for 

are obtained by using the one-band BCS formula. On the other hand, in works mm it has been shown that in a 
two-band superconductor there are two coherence lengths which are not related to the concrete bands involved in the 
formation of the superconducting ordering in a system with interband interaction: one of the lengths is diverges at 
the critical temperature ^i{T —>■ T^) —>■ 00, the second of them is a nearly constant at all temperatures ^2{T) ~ const. 
Besides it is necessary to be more accuracy at calculations of interaction between vortexes - many corrections to the 
simple GL or London theories are expected to modify the monotonically decreasing interaction potential at large 
distances, V{r) oc exp{—r/a), such that a becomes complex. This, in principle, causes an oscillating potential, whose 
first minimum may occur at large distances where the amplitude of the potential is small. Generally, as discussed in 
m, it should be taken into account dependence on the material, its purity, magnetic history, and temperature. In a 
paper m it was shown that coherence length is the same for both order parameters Ai, A2, moreover the ratio of the 
order parameters is T-independent in the GL domain, Ai(r, T)/A2(r, T) = const, with the constant depending on 
interactions responsible for superconductivity - thus the type- 1.5 superconductivity is absent. In a paper m it was 
demonstrated that close to the transition temperature, where the GL theory is applicable, the two-band problem maps 
onto an effective single-band problem with a GL parameter -|- a penetration depth -I- A^^ 

and a coherence length ^ -|- where Ki,\i, are quantities corresponding to a band i. Similar effective 

single-band GL approach also was applied in papers [nids]. The two-band GL theory has been developed in works 
[IMD where it was shown that the presence of two order parameters leads to a nonlinear temperature dependence 
of the upper and lower critical fields Hc2iT), Hci(T) and thermodynamic magnetic field Hcm(T) unlike single-band 
GL theory. In [ 32 ] the temperature dependence of the London penetration depth A(T) has been determined. These 
results are in good agreement with the experimental data for bulk MgB2 and borocarbides without any hypothesis 
about ”type- 1.5 superconductor” and ”semi-Meissner state”. 

In this paper we study two problems which, in our opinion, are important for GL theory of isotropic bulk multi-band 
superconductors: 

1 ) The coupling between the bands is represented by both the term of proximity effect Eq. @ and the term of drag 
effect - interband mixing of order parameters’ gradients: 

?7 (Vd']^V 4 ' 2 -f . ( 4 ) 

Since electron from different bands are interacting, hence, if in some a band the order parameter is spatially inhomo¬ 
geneous 'l'i(r) then in other band the order parameter must be inhomogeneous too 'L2 = 'I'2(r). If a current exists 
in one band then it drags Gooper pairs in other band. Therefore the coefficient 77 must be function of carriers’ mass 
in each band mi, m2 and the coupling e between the order parameters. As a rule the drag effect is neglected or the 
coefficient 77 is considered as an adjustable parameter. However in a work [ 24 ], where they considered Little-Parks 
effect for two-band superconductors, it has been found that the coefficient 77 is not a arbitrary quantity and a relation 
between the coefficient and effective masses of carriers exists to ensure the existence of the absolute minimum of the 
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free energy functional. In present paper we show that the drag effect plays important role in two-band superconduc¬ 
tors. Accounting of the drag effect leads to single coherence length ^ for a two-band superconductor unlike the papers 
[SHU]. Moreover the ratio of the order parameters is T-dependent Ai(r,r)/A2(r,T) = const{T), unlike the work 
m Neglecting of the drag effect leads to dynamical instability of the two-band superconductor due to violation 
of the phase relations ([^. Thus type- 1.5 superconductors are impossible. Unlike previous works we have found the 
coefficient 77 as a function of mi, m2,s. 

2 ) GL equation for a single-band superconductor (in absence of a magnetic field) is a nonlinear second-order 
differential equation. Phenomenological theory for bulk isotropic two-band superconductors has been developed in 
works [HHni, where GL equations are a set of two nonlinear second-order differential equations. Exact GL theory for 
two-band superconductors is mathematical complicated and cumbersome. Generalized set of GL equations for multi¬ 
band superconductors will be extremely complicated. Therefore approximate methods are required. In this paper we 
show that, using the result about the drag effect, the GL theory for a two-band superconductor can be reduced to the 
GL theory for an effective single-band superconductor. Generalizing this result we develop an algorithm which allows 
to reduce the free energy functional of a multi-band superconductor to the GL free energy functional of an effective 
single-band superconductor. 


II. TWO-BAND SUPERCONDUCTOR 

In presence of two-order parameters in a bulk isotropic s-wave superconductor, the GL free energy functional can 
be written as 

F = [ d^r[^ \D'i>i\^ + ^ (D+'i>+D^2 + D'i>iD+'i>+) 

J Ami 4m2 4 ^ ^ 

+ai I'I'il -l-a2|'I'2| + ^ |^i| + “^ l'f'2| + £ +'I'l^J) + ■^]) (5) 

where the differential operator are = V — (<i)o = Trhc/e is a magnetic flux quantum, H = rotA is a 

vector potential), mi^2 denotes the effective mass of carriers in the correspond band, the coefficient a is given as 
= 7i(^ ~ Tci), 7 is constant, the coefficients 61.2 are independent on temperature, the quantities e and rj describe 
interband mixing of two order parameters (proximity effect) and their gradients (drag effect), respectively. If we 
switch off the interband interaction e = rj = 0 then we will have two independent superconductors with the different 
critical temperatures Tci and Tc2 because the intraband interactions can be different gii ^ §22- There is another form 
of the coefficients they acquire constant parts consti + ji(T — T^) such that constiConst2 = and T^, is 

critical temperature of a two-band superconductor. However in this case if we switch off the interband interaction 
e = 0 consti^2 = 0, then we will have two independent superconductors with the same critical temperatures Tc. 



Figure 1: Superconductor gap parameters Ai and A 2 if the interband interaction is absent (e = 0) (dash lines) and if the 
interband interaction takes place (e ^ 0) (solid line). 

Minimization of the free energy functional with respect to the order parameters, if V'I'i_2 = 0 and A = 0 , gives 

f oi^'i -b e '^2 + ^1^? = 0 1 

02'1^2 + £'kl -|- 62 'Ll “ 0 / 


( 6 ) 
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Near critical temperature Tc we have 2 0 , hence we can find the critical temperature as a solvability condition 

of the linearized Eqs.(|6|: 


0102 — — 7172(212 — 2ci)(Tc — Tc2) — — 0. 


( 7 ) 


Solving this equation we find > Td, Tc2, moreover the solution does not depend on sign of e. The sign determines 
the phase difference of the order parameters | 4 'i|e*‘^E |4'2|e*‘^^: 


cos((/ 5 i — 1P2) = 1 if e < 0 
cos((/?i — 1 P 2 ) = —1 if e > 0 


( 8 ) 


that follows from the Eqs. and is an analogue of Eq.Q: the case e < 0 corresponds to attractive interband 
interaction ^12 > 0 , the case £ > 0 corresponds to repulsive interband interaction 512 < 0 . It should be noted that the 
interband mixing of two-order parameters £ ensures the single critical temperature of a two-band superconductor 
whilst each band has own critical temperature - Td and Tc2 if the interband interaction is absent. This fact is 
illustrated in Fig.([l|), where it is given the qualitative picture of calculations in [ 3 ]. 

Phase relations imposes restrictions on the coefficient r). Eor temperatures near and magnetic fields smaller 
than Hci, the influence of the field on modulus of the order parameters can be neglected and we assume |'hi| = 
const, | 4 ' 2 | = const. Then the wave function can be written as 4 '^ = jdtjl exp(*(^j(r)), where ‘Pj(j) are the phases of 
the order parameters. The GL free energy functional ([^ can be rewritten as 


F = 


' 27 rA 

ur[——ni ( V(/?i — 


Stoi 


h? _/ 27 rA 

- — 


4*0 
V (/32 - 


( 27 rA 

«2 V(/22 - 


8to2 

27 rA\ 


$0 

cos((pi - (^2) 


TT^ 

+£^Jnin 2 cos((pi - (P2) + 

OTT 


( 9 ) 


where ni = 2 | 4 'ip and 712 = 2|'I'2p are the densities of superconducting electrons for the corresponding bands. 
Phase relations ^ must be satisfied over the enti re volume of a superconductor: (^i(r) — :^2(r) = const, otherwise 
superconducting state will be destroyed - Eqs. (Ip). Hence the phases must change equally: 


V</ 5 i(r) = \t(p 2 {r). 


( 10 ) 


Minimizing the free energy functional @ with respect to the vector potential A we find the current J = x H: 


27 rc. f 27 rA\ 

J = 7 ,:^ h-”-1 - 

$0 4 mi \ $0 / 


f 27 rA\ 

-- 712 V(p2- 

47712 V ^0 / 




7 ?V ^1772 


2 ttK\ ( 27 rA 

- 77)+ - -w 


COs((pi - IP 2 )]- 


( 11 ) 


Let us consider a superconductor with an inner cavity. We integrate Eq. © along a closed path lying within the 
superconductor around the cavity at a distance from the cavity’s surface larger than magnetic penetration depth A. 
Hence on the path we have J = 0 and integral on the right-hand is equal to zero. Then 


+ 77-\/7li7l2 COs(v?i — ^ 32 )^ ^ VipidX + -I- 77y'7li772 COs((pi — <^2)^ W(p2dl 


ni 
mi 
27 r$ / Til 


772 


, 2rjy/riin2COs{ipi - (P 2 ) + 

4*0 XlTT-l 7712 


( 12 ) 


where $ = ^ Adi is a magnetic flux. Taking into account the functions ipi and ipi must be single-valued fVipidl = 
j> V(p2dl = 27771 , we find that the magnetic flux through the cavity takes a discrete series $ = 7 i$o like in single-band 
superconductors [El [ 24 ]. 

Let us analyze the functional The term £y/nin^cos{ipi — 1^2) < 0 always because Eq.(|^. This lowers the free 
energy. For stability of the superconducting state it is necessary that a spatial inho mog eneity of the order parameters 
enlarges the free energy. Since we have \ 7 ipi\ 7 ip 2 = (Vv^i)^ = (Vi^2)^ > 0 from Eq.( 10 l then the stability condition is 


— -I- - + 2 ri^/rv^cos(i.pi — ^p 2 ) > 0 . 

777.1 7772 


(13) 
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From the Eq. 0 we find the London penetration depth in the following form 


A-"(r) = 


47 re^ 


_|_ 2?7A/ni(T)n2(r) cos((/3i - 1 P 2 ) 


nil 


ni 2 


( 14 ) 


From this formula we can see the condition ( 13 ) ensures A^(T) > 0 when ui(T),n2(T) ^ 0 . Thus the condition ( 13 ) 
restricts the possible quantities of the parameter rj. Let Vtpi = V(/72 = 0 , that is a paramagnetic current is absent. 
Then the free energy functional takes the form 


F = 


Stt , 


+ A^ (rotH)" 


( 15 ) 


Let the field Ho is directed along the axis Oz and a superconductor are in a halfspace a; > 0 then the magnetic field 
within the superconductor are H{x) = Hoexjp (—x/X). Substituting this field in Eq.( 15 ) and integrating we have the 
free energy per unit of square: 


F = 


Stt 


( 16 ) 


We can see the smaller London penetration depth A the smaller free energy. Then from Eq. ( 14 ) it follows that such 
quantities of the parameter t], when 


r]cos{(fii — ip2) > 0 ^ Tye < 0, 


( 17 ) 


lower the free energy. 

Let us consider a case when a two-band superconductor in a normal state (0102 > e^) has contact with a metal 
in a superconducting state. Let the superconductor in a normal state occupies a halfspace x > 0 . Since in a normal 
region the order parameters are small, then minimization of the free energy functional ([^ with respect to the order 
parameters gives 




4 mi dx‘^ 




1 _l_ 

‘ ^ A 'I 


dx^ 


( 18 ) 


Eq.( 18 ) are a set of linear equations with constant coefficients. Hence we must seek a solution in a form 4 'i = 
4^2 = 1 where the quantity k has physical sense of an inverse coherence length: k = 1 /^. Then we have 


(Sl-«i)* + ('¥’)-^)v>2 = o 

(fiy, - e) lAi + (S - aj) *2 = 0 


The characteristic equation are 






mim 2 




02 


Ol 


-r] -^-r -1-2?7e -f 0102 - e = 0. 


TOi 1712 


Solutions of this equation corresponds to two coherence lengths. At T 

-2 


Tc we have 


K = 


kl = 


0-1(22 — ^ 


4 Vmi m2 ' J 

+ - 2 r]e] 

m2 ' J 

Blf 1 

4 ymim 2 ’ J 




77217722 


( 19 ) 


( 20 ) 

( 21 ) 

( 22 ) 


The first of them is fci = 0 at the critical temperature (when 0102 — = 0 ). That is the coherence length = l/|/ci| 

is diverging ^i(T —)■ T^) —)■ 00. On the contrary k2{T = T^) ^ 0 and it varies little with temperature. These length 
scales are not related to the concrete bands involved in the formation of the superconducting ordering in a system 
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with interband interaction. This result corresponds to the results in works mm obtained by microscopic approach, 
however they suggested that the intergradient interaction is absent {rj = 0 ). 

According to the method for solving of a set of linear differential equations with constant coefficients we have to 
write solutions of Eq.( 18 ) in a form 


4-2 = + C2'tp2 ^e'^^^ ’ 


( 23 ) 


where coefficients correspond to the eigenvalue ki (they must be found from Eq.( 19 ) substituting k = 

ki), the coefficients correspond to the eigenvalue ^2- Solutions ( 231 corresponds to boundary conditions 

'I'i,2(a:: —t oo) = 0 , that is fci, A:2 < 0 . A case rf' > l/mim2, when the eigenvalue ^2 is complex (the solution 4 * = 
is oscillating), will be considered below. From the first equation of Eq.( 19 ) we have: 


Ami 


ai 


For fc = fci = 0 (at T = Tc) we have 


4 '/ ^ 




( 24 ) 


( 25 ) 


Eq.( 25 l conserves the phase relations ([^: if e < 0 the condensates in different bands are in a phase, if £ > 0 the 
condensates in different bands are in antiphase (we are in a temperature region Tci,Tc2 < T < Tc hence 01,02 > 0 ). 
For k = k2 and taking into account the condition Eq. (0, which lowers the free energy of a superconductor in a 
magnetic field, we have: 


= - 


1 a2m2 + 2\rje\mim2 


rjmi 02702 + oiwi + |77e|TOim2 + \e/ri\ 


mfm2??^oi (2) 
V’l • 


( 26 ) 


In the case when the drag-effect is neglected ry = 0 we have: 


42 ) ^ ^^( 2 )^ 

emi 


We can see that Eqs .(26 2 l\ are opposite to the phase relations 
(because 77 > 0 ), for e > 


8|): when e < 0 then = const 


( 27 ) 

44 j const < 0 

0 it is analogously. This fact leads to instability of a superconducting state in a spatial 
inhomogeneous medium: any spatial inhomogeneity violates the phase relations (|^ and, consequently, suppress the 
superconducting state. For some quantities of 77 in a case rye > 0 the dynamical stability can perhaps exist, however in 
this case the London penetration depth ( 14 ) increases and, hence, the free energy ( 16 ) increases compared with a case 
rye <0. In a case rj^ > \lrfi\m2 the solution 47 = is oscillating and it does not satisfy the boundary conditions 
00) = 0 . The solution 4 ^ = could be removed supposing C2 = 0 . However the eigenvalues ki and k2 are 
derived from the intrasystem interaction and corresponds to the different length scales in the system. Consequently 
their selection by the boundary conditions is unphysical (unlike symmetric solutions k and —k one of which can be 
selected according to the boundary conditions). Thus, to ensure stability of a superconducting state and minimality 
of the free energy, the solution ki must exist only. Then from Eq.( 20 l and Eqs.( 14 |l 6 ) we can see that the coefficient 
of intergradient interaction must be 


4 = 


77117772 


Tye < 0. 


( 28 ) 


In this case we have only one eigenvalue k = ki such that k{T —>■ Tc) = 0 . 
Usind Eqs.( 21 | 2 ' 8 l and Eq.([^ we can obtain the coherence length as 




f (^ + ^+2HN) 


|aia2 — 



2lii\/s) 


( 29 ) 


At the critical temperature this coherence length diverges ^(T —>• Tc) —>■ 00 because ai{Tc)a2{Tc) — = 0 . The 
similar problem, junction between a single-band superconductor and a two-band superconductor, was considered in 
[ 26 | . where it is assumed the phase shift A 9 in the junction is zero under the condition of no current and no field. 
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However the current through the junction between two superconductors is J = JosinA0, that is the condition J = 0 
is satisfied by both Ad = 0 and A9 = tt. In [27] it has been shown the dependence of the current on the phase 
difference J oc sinA0 for the junction between a single-band superconductor and a two-band superconductor also 
takes place. Since the phase different in a two-band superconductor is either 0 or tt - Eq.([^, then proximity of a 
single-band superconductor can not change the phase relation in a two-band superconductor. 

Single coherence length allows us to represent the orders parameters in a form 'I' 2 (r) = (7(T)d>i(r), where the 
coefficient C is not function of spatial coordinates (as follows from the above, C > 0 if e < 0 and C < 0 if e > 0). 
Hence the free energy functional of a two-band superconductor ([^ can be rewritten in the form of GL functional of 
a single-band superconductor: 



where the coefficients have a form 


/ 2Tri . \ 


^ ,, ,2 B 

, ,4 

V - —A 

4- 

4-^4- -1- — 


V ^0 J 


2 

' ' Stt 


A = ai-\- 02(7^ -I- 2£C' 

B = bi+ b2C^ 


M-^ 


1 I I 2 |C| 

nil ni2 y/mim2 ’ 


(30) 


(31) 

(32) 

(33) 


and we have redesignated dt = vlti. Thus the theory of a two-band superconductor is reduced to GL theory of 
a single-band superconductor. All characteristics (coherence length, magnetic penetrations depth, GL parameter, 
critical magnetic fields, magnetization, critical currents in a wire etc.) can be found by usual GL theory. However, 
unlike GL theory, the coefficient B and the effective mass M are functions of temperature since the coefficient G is a 
function of temperature. 

Now we should find the coefficient C. Let us substitute d >2 = GTi in Eq.([^: 


f di eC -|- — 0 1 

\a2C + e + b2C^'i>l = 0 j 


(34) 


If r —>■ Tc then the equations can be linearized. In this case we have solutions C = —aije or C 
critical temperature we can use Eq.Q, that is |e| = 1 / 0102 . Then the solution becomes unique: 



£ < 0 
£ > 0 


—£/a 2 . Near the 


(35) 


This approximation expresses the fact that relation between the order parameters is determined by the single-band 
critical temperatures Tci,T’c 2 ' if Td > then > A 2 - Fig. 0- 

Using Eqs.( 31|32|33|34 ) we can find main characteristics of a superconductor as in the usual GL theory. A coherence 
length: 




^2 

4M|A| 



4 2 |di«2 - kl 


i/mim 2 


a magnetic penetrations depth: 


a GL parameter: 


Mc^B 

87re2|A| 




bl^ 
^ ai 


a-2 




yjm-im2 


n = 


A 




2 

y/mim2 


) 


(36) 


(37) 


(38) 


We can see that the GL parameter is a function of temperature unlike single-band GL theory. However this dependence 
is type that is little varying function of temperature if T ^ Tci,Tc 2 - It should be noticed that this approximation 
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is correct if T > Tci,Tc 2 only. We can extrapolate the obtained expressions for all temperature. To do it we can 
suppose Oi = 7 i(T'c — Td) = const, then M = const, B = const, however it is necessary to expand the expression 
— jel in powers of T — T^: 


yja^- |e| 


7i 72 — Tel ~ Tdl) 

2^01 a2 


(r-T,) 


7i72 

2V“1«2 



7i72 {2Tc - Tel - 
40102 



(T -T,f + ... 


(39) 


Thus in the functional for a two-band superconductor the coefficient ^ is a power series of (T — Tc) unlike the GL 
functional for a single-band superconductor. From this fact a nonlinear temperature dependence of the upper critical 
field follows (hear fi, /2 are some coefficients): 


<i>o (Tc — T) single-band GL theory 

27r^2 fi {Tc — T) + f 2 {Tc — T)^ -|- ... two-band GL theory 


(40) 


that is consistent with experimental data (in bulk LuNi 2 B 2 C, MgB 2 ) in pHl - IHH] and theoretical results in [TTIIT!^ . where 
it has been shown that the presence of two order parameters for two bands yields a nonlinear temperature dependence 
of Hc 2 {T) in the vicinity of the critical temperature unlike the single-band s-wave BCS theory and GL theory. It 
should be noted that this difference can be a cause of strong enhancement of Hc 2 {T) (up to ten-fold increase) in dirty 
two-gap superconductors, that, as noted in |31] . is result from the anomalous upward curvature of Hc 2 {T). For the 
lower critical field Hd and the thermodynamic magnetic field Hem we have analogous expansion because 


Hr. = 


^>0 


27rA^ 
4>o 


In K cx l^l — \Jaia 2 




(41) 

(42) 


that demonstrates nonlinear temperature dependence too and correlates with theoretical results of nanni. Let the 
carriers have different effective masses in different bands, for example mi ^ m 2 . From Eq.(|33| we c an see that 
the two-band effective mass M is determined mainly by the smaller mass m 2 . From Eqs. ( 36|37|40|41 ) we can see 
the critical fields Hci and Hc 2 depend on the effective mass as Hci oc 1/M, Hc 2 oc M. Hence, as noted in [TU1121) . 
the critical fields are determined mainly by the smaller mass m 2 , while the contribution from the lager mass can be 
neglected. 


III. MULTI-BAND SUPERCONDUCTOR 

Using results of previous section we can generalize the above-described method for two-band superconductors to 
multi-band superconductors. In presence of n order parameters in an isotropic s-wave superconductor, the free energy 
functional can be written as 


F= f |UvI-,|2 + a, + I 14/,|4 

n ^2 

j=2,j>i 

Critical temperature can be found from an equation: 


ai 

^12 ■ 

■ ^ In 


^12 

02 . 

■ E2n 

= 0 , 


S2n ■ 

O'n 



(43) 


( 44 ) 
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which is analog of Eq.Q. However we should notice that in general case the symmetry e o —e for critical temperature, 
like in the two-band case, is absent. If all Sij < 0 but some > 0, suppression of superconductivity is possible. We 
will consider a case of attractive interband interaction only, that is all Eij < 0 . 

Following our scheme we should find coefficients of the intergradients interaction Tyij and the coherence length 
Equation for the coherence length = 1/fc^ is 


- ai ^r]i2k‘^ - £12 • ■ • - £in 

- £12 4^^^ - ^2 • ■ • “ £2n 

^Vlnk'^ - £ln ^V2nk'^ “ £2n • ■ • 4 ^^^ “ «n 

= fnimi,Vij)k‘^"" + -b ... -b /i(mj,?7y )fc^ + /g = 0, (45) 


which is analog of Eq.(20). At T = Tc we have /o = 0. The coefficients rjij must be such that the functions 
fn = f{n-i) = • ■ • = /2 = 0 , then the coherence length is 


h 


By analogy of (28 1 and using the condition Eu < 0 we can suppose 


1 


= 


^yrn^rrij 


(46) 


(47) 


In the next step we should to represent the orders parameters in a form il '2 = < 112 (T)!?!, 413 = C' 3 (T)'I'i,..., = 

Cn(T)4'i. Then the free energy functional of a multi-band superconductor (43) takes the form of the GL functional 
(30) of a single-band superconductor with coefficients 


A — ai + ciiCf + 2 £iiC'i -b 2 EijCiCj 

i—2 i—2 i—2 j—3,j>i 

n 

B = h+J2hCf 

2=2 


M-^ 


1 

mi 


+ E 


2 = 2 


c? 


rrii 


n n n 

-b 2 ^ r]iiCi + 2 ^ ^ Vij ^2 Cj 

i—2 i—2j—3J>i 


(48) 

(49) 

(50) 


Linearized equations for ( 72 , ( 73 ,..., C„ are 


Ol -b £ 12^2 -b ...EinCn — 0 
£12 + 02(72 -b ■..E 2 nCn = 0 

£lra + £23^2 -b ...a„(7„ = 0 


(51) 


which have to be solved taking into account Eq.(44) so that the solutions are unequivocal (as we have shown in the 


two-band case). However we can use an approximate method. In the two-band problem we supposed the coefficient 
C = 


^ for 4(2 = (74'i, that is relation between the order parameters is determined by the single¬ 

band critical temperatures Td, Tc2 - if Tci > Tc 2 then Ai > A 2 . This fact can be used for the coefficients Ci in the 
multi-band problem, where we can suppose: 



(52) 


Substituting Eq.(52) in Eq.(48) and reducing to a common denominator we have 

nydT 


A = 


n :=2 


/ai 


-f {ai^ Eij ), 


(53) 
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where 


n 2 2 /- '*■ 

f "h ^ ^ x/“t“ ^ ^ ^ ^ x /: (^4) 

i=l i=2 fc=2,fc/i i=2 j=3,j>i k=2,k^i,k^j 

The critical temperature is such a temperature when f(T = Tc) = 0. As in the two-band problem we can extrapolate 
the obtained expressions for all temperature. To do this we have to suppose = 7 i(T'c — Td) = const, then 
M = const, B = const, however it is necessary to expand the expression f{ai,eij) in powers of T — Tc- Thus the 
multi-band problem is reduced to the single-band problem with the effective mass M, however the coefficient A is 
power series of (T — T^) unlike the GL free energy functional. 


IV. RESULTS 


In this work we have shown that the term of the drag effect rj (V'kf Vd >2 + V'kiV^'J) in the free energy functional 
of an isotropic bulk two-band superconductor plays important role and the restrictions for the coefficient 77 exist. 


If the coefficient is 77 ^ = 


and it’s sign is opposite to the sign of the coefficient in the term of the proximity 


^ '^2 + that is rje < 0 , then this leads to a single coherence length which diverges at the critical 


effect £ ('kf 

temperature ^(T —>■ T^) —> 00 , and a single GL parameter. This quantity ensures the stability of a superconductor 
state and the least possible free energy in this case. Other quantities of the coefficient or neglecting of the drag effect 
77 = 0 leads, at first, to the existence of two coherence lengths, where one of them diverges at the critical temperature 
while the second length is finite at all temperatures. Secondly, it leads to the dynamical instability (suppressing of 
a superconducting state if the order parameters are spatial inhomogeneous) due to violation of the phase relations 
(^). These results mean that the isotropic bulk type-1.5 superconductors are impossible. 

It should be noticed that these results are obtained in the GL domain only. Hence it can be supposed that at low 
temperatures the disproportion 'k 2 (i",?") ^ C'(T)'I'i(r,T) can takes place, that is there are two different coherence 
lengths 7 ^ ^ 2 - However this fact means that the order parameters have different gradients V4'i(r) 7 ^ V'I' 2 (r). Since 
the order parameters are |'I'i|e*‘^L |^' 2 |e*‘^’^, then the different gradients can lead to violation of the equality (lOl, 
hence to violation of the phase relations (|^. Thus the state with different coherence lengths is dynamically unstable. 

The approximate method for solving of set of GL equations for an isotropic bulk multi-band superconductor has 
been developed. Using the results about the drag effect we have shown that the free energy functional for a two- 
band superconductor can be reduced to the GL functional for an effective single-band superconductor. This effective 
superconductor is characterized with some an effective mass of carriers (as a function of 777,1,7712,77) and a coefficient 
at I'kp as a power series of (T — T^) in the vicinity of the critical temperature. This temperature dependence causes 
nonlinear dependence of upper and lower critical fields Hc 2 , Hd, thermodynamical magnetic fields Hem on temperature 
unlike the single-band GL theory. Generalizing this result we have developed an algorithm which allows to reduce 
the free energy functional of a multi-band superconductor to the effective GL free energy functional of a single-band 
superconductor provided that all interband interactions are attractive. 


Acknowledgments 

The work is supported by the project #0112U000056 of the National Academy of Sciences of Ukraine. 


[1] N. Kristoffel, P. Konsin and T. Ord, Rivista del nuovo cimento 17, 1 (1994) 

[2] Oleg V. Dolgov, Igor I. Mazin, David Parker, Phys.Rev.B. 79 060502(R) (2009) 

[3] G. Litak, T. Ord, K. Rago and A. Vargunin, Asta Physica Polonica A 121, 747 (2012). 

[4] A. Ciechan, Karol I. Wysokiski, Phys.Rev.B. 80 224523 (2009) 

[5] A. J. Leggett, Prog. Theor. Phys. 36, 901 (1966) 

[6] S.G. Sharapov, V.P. Gusinin, H. Beck, Enr.Phys.J.B 30, 45 (2002) 

[7] A. Gnrevich, V. M. Vinokur, Phys. Rev. Lett. 30, 047004 (2003) 

[8] V. Moshchalkov, M. Menghini, T. Nishio, Q. H. Chen, A. V. Silhanek, V. H. Dao, L. F. Chibotarn, N. D. Zhigadlo, and 
J. Karpinski, Phys. Rev. Lett. 102, 117001 (2009) 

[9] T. Nishio, V. H. Dao, Q. Chen, L. F. Chibotarn, K. Kadowaki, and V. Moshchalkov, Phys. Rev. B 81, 020506(R) (2010) 

[10] E. Babaev, J.M. Speight, Phys. Rev. B 72, 180502 (2005) 

[11] E. Babaev, J. Carlstrom, J.M. Speight, Phys. Rev. Lett. 105, 067003 (2010) 







11 


[12] E. J. Carlstrom, E. Babaev, M. Speight, Phys. Rev. B 83, 174509 (2011) 

[13] E.H. Brandt, M. Das, Journ. Supercond. Novel Magnetism, 24, 57 (2011) 

[14] T 6rd, K. Rago and A. Vargunin., J. Supercond. Nov. Magn. 25, 1351 (2012) 

[15] V.G. Kogan, J. Schmalian, Phys.Rev.B 83 054515 (2011) 

[16] Jani Geyer, Rafael M. Fernandes, V.G. Kogan, and lord Schmalian, Phys. Rev. B 82, 104521 (2010) 

[17] I.N. Askerzade, A. Gencer, N. Giiglii, Supercond. Sci. Technol. 15, L13 (2002) 

[18] I.N. Askerzade, A. Gencer, N. Giiglii and A. Kilig, Supercond. Sci. Technol. 15, L17 (2002) 

[19] I.N. Askerzade, Physica C 397 99 (2003) 

[20] I.N. Askerzade, Physics-Uspekhi 49 1003 (2006) 

[21] Iman Askerzade, Unconventional superconductors. Anisotropy and multiband effects., Springer, 2012, p.28 

[22] I.N. Askerzade, A. Gencer, Solid State Communications 123, 63 (2002) 

[23] H. Doh, M. Sigrist, B.K. Cho, S.I. Lee, Phys. Rev. Lett. 83, 5350 (1999) 

[24] Yu.S. Erin, S.V. Kuplevakhskii and A.N. Omelyanchuk, Low Temp. Phys. 34, 891 (2008) 

[25] M.E. Zhitomirsky, V.H. Dao, Phys.Rev.B 69, 054508 (2004) 

[26] Y. Tanaka, Journal of the Physical Society of Japan 70 No.10, 2844 (2001) 

[27] A. Brinkman, A. A. Golubov, M. Yu. Kupriyanov, Phys.Rev. B 69 214407 (2004) 

[28] J. Freudenberger, S.-L. Drechsler, G. Fuchs, A. Kreyssig, K. Nenkov, S.V. Shulga, K.-H. Muller, L. Schultz, Physica C 
306, 1 (1998) 

[29] K.-H Muller, G. Fuchs, A. Handstein, K. Nenkov, V.N Narozhnyi, D. Eckert, J. Alloys Compd. 322 LIO (2001) 

[30] S.L. Budko, C.Petrovic, G.Lapertot, C.E. Cunningham, P.C. Canfield, M-H. Jung, and A.H. Lacerda, Phys. Rev. B 63 
220503(R) (2001) 

[31] A. Gurevich, Physica C 456, 160 (2007) 



